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We consider measurements, described by a positive-operator-valued measure 
(POVM), whose outcome probabilities determine an arbitrary pure state of a D- 
dimensional quantum system. We call such a measurement a pure-state information- 
ally complete (PS I-complete) POVM. We show that a measurement with 2D — 1 
outcomes cannot be PS I-complete, and then we construct a POVM with 2D out- 
comes that suffices, thus showing that a minimal PS I-complete POVM has 2D out- 
comes. We also consider PS I-complete POVMs that have only rank-one POVM 
elements and construct an example with 3D — 2 outcomes, which is a generalization 
CN ■ of the tetrahedral measurement for a qubit. The question of the minimal number of 

elements in a rank-one PS I-complete POVM is left open. 

> 

m ■ I. INTRODUCTION 

An important technical requirement for developing coherent quantum technologies is the 
ability to assess how well one can prepare or create a particular quantum state. This assess- 
ment is carried out by making suitable measurements on a sequence of identically prepared 
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quantum systems. The measurements are chosen so that the probabilities of the mea- 
surement outcomes suffice to determine the state; the probabilities are estimated from the 



outcome frequencies observed in the measurements. A set of measurements whose outcome 
probabilities are sufficient to determine an arbitrary quantum state is called informationally 
complete 0, . The process of determining the quantum state from the results of a sequence 
of these measurements is called quantum state tomography. 

Given a set of measurements that are informationally complete, they can be amalgamated 
into a single generalized measurement, consisting of a coin flip to select a measurement 
from the set, followed by an application of that measurement. The statistics of such a 
generalized measurement are described by a positive-operator-valued measure (POVM) 0, 
which consists of positive operators E c that are labeled by an index c that runs over the 
possible outcomes. These operators satisfy a completeness condition, 

c=l 

The operators E c , satisfying < E c < I, are called POVM elements. The probability for 
outcome c, given system state p, is 

P c = tr(pE c ) . (2) 

The completeness condition and the unit-trace normalization of p guarantee that the prob- 
abilities are normalized to unity. 

An informationally complete set of measurements can thus always be reformulated as 
a single informationally complete POVM (IC-POVM), i.e., a POVM having the property 
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that the relations giving p c in terms of p can be inverted to give p in terms of the outcome 
probabilities. A normalized density operator p is specified by D 2 — 1 real numbers, since 
an arbitrary Hermitian operator on a D-dimensional quantum system is specified by D 2 
independent real numbers, and this number is reduced by one for a normalized density 
operator because it has unit trace. Since the relations (j2J are linear, it is a simple matter 
of linear algebra to conclude that an IC-POVM must provide D 2 — 1 independent relations. 
This might lead one to think that an IC-POVM can get by with D 2 — 1 elements, but the 
normalization of the outcome probabilities means that one of the relations (J2J) is redundant, 
so an IC-POVM must have at least D 2 elements. 

An equivalent, often more useful way of thinking is to disregard the normalization of the 
density operator. Then, like any Hermitian operator, p is specified by D 2 real numbers, the 
additional number being tr(p). By virtue of the completeness condition, any POVM provides 
tr(p) as the sum of the (unnormalized) probabilities This way of thinking bypasses the 
waffling about normalization conditions and goes directly to the point: an IC-POVM must 
provide D 2 linearly independent relations and thus must have at least D 2 POVM elements. 

Informational completeness for arbitrary quantum states is thus equivalent to the ability 
to reconstruct any Hermitian operator H from the operator inner products tr(HE c ). This 
means that the problem of finding informationally complete POVMs is reduced to finding 
positive operators F c that span the vector space of Hermitian operators. If these operators 
don't satisfy the completeness condition (0), one replaces them with the operators E c = 
G-^-F'G- 1 / 2 , where G = is trivially a nonsingular operator. The new operators are 

positive, satisfy the completeness condition, and are obviously informationally complete. 
Using this procedure, it is easy to construct minimal IC-POVMs, i.e., POVMs having D 2 
elements, from sets of D 2 linearly independent, positive operators. Indeed, it is easy to 
construct minimal IC-POVMs for which all the POVM elements are rank-one, i.e., multiples 
of projectors onto pure states 0. Moreover, there is evidence, both analytical and numerical, 
that in all finite dimensions there are minimal, rank-one IC-POVMs, called symmetric, 
informationally complete POVMs (SIC-POVMs), for which the operator inner products of 
all pairs of POVM elements are the same Recast as ensembles of states, SIC-POVMs 
have applications to quantum key distribution 0, and they have been shown 0] to be the 
ensembles that are most quantum according to a measure of quantumness introduced by 
Fuchs and Sasaki 0. 

In this paper we consider the pure-state version of informational completeness, i.e., re- 
construction of an arbitrary pure state p = \ip){ip\ from the outcome probabilities 

p c =^\E c \^). (3) 

We formalize the notion of pure-state information completeness in the following definition. 

Definition (PS I-completeness). A pure-state informationally complete (PS I- complete) 
POVM on a finite-dimensional quantum system is a POVM whose outcome probabilities are 
sufficient to determine any pure state (up to a global phase), except for a set of pure states 
that is dense only on a set of measure zero. 

The intent of the last clause is to require the inversion to be unique for generic pure states. 
This clause says that any pure state outside a set of measure zero is surrounded by an open 
ball of states all of which are uniquely determined by the outcome probabilities. 

Though some of the thinking and some of the mathematical techniques developed for IC- 
POVMs can be transferred to a study of PS I-complete POVMs, there is a critical difference, 
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which makes PS I-completeness more difficult to analyze: the relation between the outcome 
probabilities and a density operator is linear, whereas the relation between outcome prob- 
abilities and pure states is quadratic. The problem of IC-POVMs thus lies squarely within 
linear algebra, whereas PS I-completeness must be tackled using other tools. 

The first point to make about PS I-completeness arises from a simple counting argument. 
A pure state is specified by D complex amplitudes, corresponding to 2D real numbers, 
but the number of independent real numbers is reduced to 2D — 2 by the normalization 
condition, (ip\ip) = 1, and by the fact that a pure state can be multiplied by a global phase 
without changing any of its predictions for probabilities. Because of the redundancy due 
to the normalization of the POVM probabilities, a PS I-complete POVM must contain at 
least one more element than the 2D — 2 real numbers to be determined. This suggests 
that a PS I-complete POVM might get by with only 2D — 1 elements. This argument 
can be profitably rephrased in terms of unnormalized pure states, which are specified by 
2D — 1 real numbers, giving immediately that a PS I-complete POVM must have at least 
2D — 1 elements. Unlike the similar counting argument for density operators, however, 
the counting argument for PS I-complete POVMs is only suggestive, because linear algebra 
can't be applied to give rigorous conclusions. The counting argument leaves us only with 
a prejudice that the minimal number of elements in a PS I-complete POVM ought to be 
2D — 1 or a little more. 

The point of this paper is to provide some rigor. We establish that 2D, not 2D — 1, is the 
minimal number of elements in a PS I-complete POVM. To show this, we first prove that a 
POVM with fewer than 2D elements — in particular, one with 2D — 1 elements — cannot be 
PS I-complete (Sec. Illl) ). and we then construct an example of a PS I-complete POVM with 
2D elements (Sec.UVjl. 

Asher Peres considered a particular version of this problem in his classic 1993 quantum- 
mechanics textbook |3j . He noted that the probabilities in two complementary bases should 
be sufficient to determine a pure state up to a finite set of ambiguities. If we let {\ej), j = 
1, . . . , D}, denote an orthonormal basis, the complementary basis consists of the vectors 



fk) = ^= y £\e i )^ D ■ (4) 



D 

D 1 1 

These complementary bases are discrete analogues of the position and momentum bases of a 
particle moving in one spatial dimension. A consequence of our work is that measurements in 
two complementary bases cannot be PS I-complete. To see this, notice that the amalgamated 
POVM for the two measurements has 2D POVM elements, 

771C _ j 2 l e c)( e c| > C — 1, . . . , D, . . 

~ UI/ c -dX/c-d| , c = D + l,...,2D. {0> 

The factor of 1/2 takes care of the coin flip that chooses with equal probability between 
measurements in the two bases. The probabilities associated with the last POVM element 
in each basis, p D = §| | 2 and p 2D = ||(/d|^)| 2 , provide no information because of the 
separate normalization of the measurements in the two bases. Without losing any informa- 
tion, we can combine these two POVM elements into a single element, ^(\en) {eD\ + \fD) (fo\) , 
thus reducing the number of POVM elements to 2D — 1. Our proof in Sec. II I II shows this 
cannot be PS I-complete; attempts to find a pure state from the outcome probabilities must 
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at least suffer from the finite ambiguities mentioned by Peres. Following Peres, we illustrate 
these ambiguities in two dimensions in Sec. ITT1 

The problem considered by Peres is often called the Pauli problem, as it is the finite- 
dimensional version of a question posed by Pauli in a footnote to his quantum-mechanics 
article in Handbuch der Physik jjj: can the wave function ip(x) of a particle be determined 
(up to a global phase) from the position and momentum probability distributions, |^(x)| 2 
and \ip(p)\ 2 7 The answer to Pauli's question is a definitive no: there are many wave func- 
tions not determined by their position and momentum distributions and thus said to be 
Pauli nonunique. Much work has been devoted to investigating Pauli nonunique states (see 



Refs. |2|, |lOj, |TJj and references cited therein), but even now there does not seem to be a com- 



plete characterization of such wave functions. A simple example of a Pauli nonunique state is 
an eigenstate of parity, ij)(x) = ±if)(—x) (and ip(p) — ±4>(—p)), that is not time- reversal in- 
variant, i.e., <j)(x) = ip*(x) ^ e ia ip(x); such a wave function satisfies = 4>{~P) = =t<Kp)) 
so ip{x) and <p(x), though distinct states, have the same position and momentum distribu- 
tions. Discussion of work on Pauli nonunique states lies outside the scope of our paper, 
but we do note the interesting result of Corbett and Hurst 0| that Pauli nonunique wave 
functions are dense in the space of wave functions, a result that motivates the denseness 
restriction in our definition of PS I-completeness. 

Previous work on PS I-completeness for finite-dimensional quantum systems k| 3 
has considered the system to be a spin-s particle (D = 2s + 1) and has focused mainly 
on reconstructing a pure state from the probabilities for spin components along several 
directions. The angular-momentum algebra and the rotation group play essential roles in 
these considerations. We review some of these results during the course of our discussion and 
point out how they are related to our work and, in particular, how they can be rephrased 
in terms of a single POVM. 

Our paper is organized as follows. In Sec. m we introduce informationally complete and 
PS I-complete POVMs in two dimensions and note that the two-dimensional case provides 
little guidance for generalizing PS I-completeness to higher dimensions. In Sec. Illl| we 
prove that the outcome probabilites of a POVM with 2D — 1 elements are insufficient to 
determine a general pure state in D dimensions, leading to, at best, a two-state ambiguity. 
In Sec. IIVI we construct a PS I-complete POVM with 2D elements, and in Sec. we turn 
to the question of rank-one PS I-complete POVMs and construct one with 3D — 2 POVM 
elements. Section fVll gives a brief summary of open questions, while Sec. IVHI closes with a 
poetic description of one of these questions, dedicated to Asher Peres on the occasion of his 
70th birthday [11. 



II. INFORMATIONALLY COMPLETE POVMS IN TWO DIMENSIONS 

We can gain some insight into IC-POVMs and PS I-complete POVMs by looking at the 
case of a two-dimensional quantum system (qubit), where the Bloch representation of states 
and operators permits us give a complete characterization of information completeness. 

The Bloch representation of an arbitrary POVM element for a qubit is 

E = al + bn ■ <x , (6) 

where n is a unit vector in M 3 and a and b are nonnegative real numbers satisfying b < a 
and b < 1 — a, to ensure that E and I — E are positive operators. Rank-one POVM elements 
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have a = b. A POVM is made up of n such elements, 

E c = a c l + b c n c -a , c=l,...,n, (7) 
satisfying the completeness condition |Tj) . which now becomes two conditions, 

1 = $> C , (8) 

c 

= ^6 c n c . (9) 

c 

A POVM is informationally complete if and only if the POVM elements E° span the 
four- dimensional space of Hermitian operators. Translated to the Bloch representation, this 
says that a POVM is informationally complete if and only if the vectors b c n c span R 3 . A 
direct way to see this in the Bloch representation is to note that for an arbitrary density 
operator p = ~(I + P ■ a), with polarization vector P (|P| < 1), the outcome probabilities 
are given by 

p c = ti{pE c ) =a c + b c n c ■ P ; (10) 

to reconstruct an arbitrary polarization vector, the vectors b c n c must span M 3 . 

The minimal number of elements in an IC-POVM is four; one sees this directly in the 
Bloch representation by noting that the vectors b c n c must span M 3 and also satisfy the 
condition (jH)), implying that there must be at least four vectors. For a minimal IC-POVM, 
it is also easy to see from condition Q that any three of the four vectors b c n c must be 
linearly independent and thus span R 3 (and this means that none of the b c can be zero). 

An example of a minimal, rank-one IC-POVM, which we make use of in Sec. is given 
by the tetrahedral measurement, which is specified by the four unit vectors 



n 1 = (0,0,1), 




with a c = b c = 1/4, c = 1,2,3,4. This POVM is also a SIC-POVM by virtue of the 
symmetric placement of the Bloch vectors, which connect the origin to the vertices of a 
tetrahedron whose apex is at the north pole of the Bloch sphere. 

A PS I-complete measurement is one such that the outcome probabilities (fTUj) determine 
uniquely a generic pure state \m)(m\ = |(I + m ■ cr). It is immediately clear that to re- 
construct an arbitrary unit vector m, the vectors b c n c must span IR 3 , so a PS I-complete 
measurement in two dimensions is always informationally complete for all states, pure or 
mixed. For this reason PS I-completeness for qubits provides little guidance for what hap- 
pens in higher dimensions: in two dimensions, the 2D POVM elements required for PS I- 
completeness provide the D 2 elements necessary for informational completeness for all states, 
whereas in higher dimensions, there is a yawning gap between 2D and D 2 . 
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It is worth stressing the way a three-element POVM fails to be PS I-complete. The 
condition (JUJ) implies that the three vectors b c n c are linearly dependent and thus span at 
most a plane. The outcome probabilities (fTU|) provide no information about the component of 
m orthogonal to the plane. When the vectors b c n c do span a plane, the outcome probabilities 
determine the projection of m onto this plane, which specifies everything about m except 
the sign of the component orthogonal to the plane. This two-fold ambiguity is the source of 
the POVM's incompleteness for pure states. 

A symmetric example of a three-outcome, rank-one measurement is the trine measure- 
ment on a qubit 15], given by the unit vectors 

n 1 = (1,0,0), 

with a c = b c = 1/3, c = 1,2,3. These vectors point from the origin to the vertices of 
an equilateral triangle lying in the equatorial plane of the Bloch sphere. The outcome 
probabilities for the trine determine the equatorial component of the Bloch vector for a 
pure state, but can't resolve whether the Bloch vector is in the northern hemisphere or the 
southern hemisphere, thus leaving a two-state ambiguity. 

III. NECESSITY OF 2D OR MORE OUTCOMES 

In this section we show that the outcome probabilities of a POVM with 2D — 1 outcomes 
are insufficient to determine a generic pure state. In particular, we prove the following 
theorem. 

Theorem. The outcome probabilities p c = (tp\E c \ilj) of a POVM with 2D — 1 measurement 
outcomes are insufficient to determine a generic pure state in D dimensions. 

We prove the theorem by showing that any allowed distribution of outcome probabilities, 
except a set of distributions of measure zero, is consistent with at least two possible pure 
states. The proof occupies the remainder of this section. 

We begin by picking a particular POVM element, say E 1 , and writing it in its orthonormal 
eigenbasis {\j), j = 0, . . . , D - 1}: 

D-l 

E l = Y,\\3)- (13) 

3=0 

We order the eigenvalues from smallest to largest, i.e., < Ao < Ai < . . . < \d-i- We can 
assume that \d-i > Ao, because otherwise E l would be a multiple of the identity operator 
and would provide no useful information beyond that always contained in the normalization 
constraint. Any pure state can be expanded in this eigenbasis as 

D-l 

M=r |0> + 5>|j>. (14) 
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We use the global phase freedom to make r real, but we do not require to be normalized. 
Instead we consider the norm to be an extra parameter of the state, which is determined by 
the sum of the (unnormalized) outcome probabilities ©, 

2D-1 

E p c = (V#> = R 2 , (is) 

c=l 

by virtue of the completeness condition (P) for our POVM. 

The amplitudes Cj are arbitrary complex numbers, which we now write in terms of their 
real and imaginary parts, i.e., Cj = Xj + iyj. The variables necessary to define the state are 
now the 2D — 1 real co-ordinates r , £i, . . . , x^d-Ii Vi, ■ ■ ■ , V2D-1-, each of which is free to 
take on any real value. We find it convenient to put all these co-ordinates into a single real 
vector £ 6 IR 2 " 0-1 , whose components are 

£o = r , 

= x i and €j+D-i = Vj , j = 1, • • • , D - 1. (16) 

Notice that an inversion through the origin, i.e., £ — > produces a global sign change in 
When r = £o 7^ 0, we could remove the global phase freedom entirely by making £ 
positive, but we choose to retain this two-fold ambiguity so that the space we are dealing 
with is the entirety of IR 2Z)_1 . 

Each of the 2D — 1 outcome probabilities can now be written as a quadratic form 

P c =^\E c \^) = J2CrM^^ sl (17) 

r,s 

where the quantities M^ s are the matrix elements of a real, symmetric, positive (semi- 
definite) matrix M c . A surface of constant p c > is a 2{D — l)-dimensional (possibly 
degenerate) hyperellipsoid S p c, which extends to ±oo in a number of dimensions given by 
2D — 1— rank(M c ) and has a hyperellipsoidal cross-section in the remaining dimensions, num- 
bering rank(M c ). Notice that by construction the first outcome probability has a diagonal 
matrix M 1 , 

D-1 D-1 

P 1 = Aor 2 + A il c ;| 2 = Ao£o + E + %+D-i) , (18) 

3=1 3=1 

and is thus a hyperellipsoid aligned with the co-ordinate axes. 

We now replace two of the outcome probabilities by equivalent constraints. First we 
replace the last outcome probability (c = 2D — 1) by the normalization condition (fTB^) . 

2D-2 

R 2 = (V#> = E £ > (19) 

r=0 

something we can do because p 2D ~ l = R 2 — Y^ 2 =\ 2 P c - A surface of constant norm R 2 is, of 
course, a sphere Sr of radius R in IR 2jD_1 . Each physical state with £ ^ is represented 
twice on this sphere of radius R, once in the northern hemisphere (£o > 0) and again, with 
its sign reversed, at the antipode in the southern hemisphere (£o < 0). States on the equator 
(£o = 0) retain the full global phase freedom. 
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FIG. 1: Transversal and nontransversal intersections. The two ellipses shown in each sub figure are 
submanifolds of M. 2 . a) shows a nontransversal intersection of two ellipses, and b) is a close-up of 
one of the intersection points. The intersection in b) is nontransversal because the tangent vectors 
at the point shown are the same; they therefore span only a one-dimensional subspace of M. 2 . c) 
shows a slight perturbation of a), and by looking at the close-up of the intersection in d), we see 
that the intersection is now transversal. This is because at all points of intersection, the tangent 
vectors of the two submanifolds span M 2 . 

The second modification is to replace the first outcome probability (fTKJ) with an equivalent 
constraint in which £o does not appear: 

D-l 

p 1 - X R 2 = - Ao)(£| + ^-i) = P > • (20) 

That P > follows because Ao is the smallest eigenvalue. Moreover, we are guaranteed that 
P is not always zero because Xd-i > Ao- All this means that a surface of constant P > is 
a hyperellipsoid Sp, aligned with the £ r axes, which runs off to ±oo along the £o axis (and 
perhaps other axes) and which has an ellipsoidal boundary in one or more co-ordinates £ r . 

Now consider a particular set of values for R 2 , p 2D ~ 2 , . . . ,p 2 , and P, and imagine applying 
the constraints imposed by these 2D — 1 measured values in the order listed. The first step 
in the process yields the submanifold X\ = S R 2 = Sr. The set X m obtained after m steps is 
the intersection of I m -i with the hyperellipsoidal submanifold S p 2D- m , 

The last step intersects X2D-2 with Sp to yield the set X2D-1 = I2D-2 H £p. We note that 
for all m, X m is a compact set because the process starts with the sphere Sr, which is itself 
compact. For PS I-completeness the final set I2D-1 must consist generically of just the two 
antipodal points representing a particular pure state. 

Having set up the entire problem, we now need a bit of differential geometry about the 
intersection of submanifolds. Given a manifold R of dimension dp, two submanifolds, M 
and N, of dimensions du and dpj, intersect transversally if at each point of intersection 
V, the vectors in the tangent spaces to M and N, Tm(V) and T^(V), span the entire 
tangent space to R at V (see Fig. |TJ). A transversal intersection M fl is a submanifold 
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FIG. 2: The cylinder is the hyper ellipsoidal submanifold £p with many dimensions suppressed. 
Recall that £p extends to infinity, though it is drawn as finite here. The union of the closed curves 
A, B, C, and D is the submanifold T.2D-2- The convention for dashed and dotted lines is that 
lines that are behind one layer of the cylinder are dashed, and lines that are behind two layers are 
dotted. The points where I2D-2 intersects £p (denoted as black dots) are the (zero-dimensional) 
submanifold Tid-\- Notice that the entire diagram is invariant under inversion through the origin. 
The shaded disk Xq, which is the part of the £0 = hyperplane that is in or on the interior of 
the £p submanifold, also has inversion-symmetric points of intersection with Z2D-2- The part of 
£p with £0 > 0, denoted in the text as £ P + \ has multiple intersections with Z2D-2, implying that 
Eqs. lfT7|) cannot be inverted to reconstruct a unique pure state \tp). 

of dimension dMnN = &m + d,N — dp, since this relation holds for the dimensions of the 
tangent spaces at each poin t of intersection. A fundamental theorem of differential geometry, 
called Sard's theorem [l6|, asserts that if two submanifolds with a nonempty intersection 
intersect nontransversaliy, they can be perturbed slightly to intersect transversally. Figure Q 
illustrates the difference between transversal and nontransversal intersections and how a 
perturbation of a nontransversal intersection yields a transversal intersection. 

We now apply these ideas to our situation. Notice that if the successive intersections 
in our sequence are transversal, then each intersection produces a new submanifold with 
dimension decreased by one, dj m = dx m _ 1 + (2D — 2) — (2D — 1) = dj m _ 1 — 1, so the 
dimension of I m is 2D — 1 — m. We now argue that if the measurement is to be PS I- 
complete, a sequence of such transversal intersections must occur generically, i.e., for all 
outcome probabilities p c except for a set of measure zero. 

Suppose instead that the intersections are transversal up to the mth step, at which step 
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intersect nontransversally. There are two ways this could happen. First, 
it could be that part or all of T m _i lies in S p 2D- m , so that part or all of the intersection fails 
to have dimension reduced by one. If this situation were generic, it would mean that the 
final intersection would be greater than zero-dimensional, thus not specifying points, so the 
POVM could not be PS I-complete. Second, it could be that part or all of the intersection 
has dimensionality reduced by more than one, as in the case of two ellipsoids touching at 
tangent points. This is the situation depicted in Figure EJaj if the ellipses were replaced 
with surfaces of revolution about the vertical axis in the diagram. If this situation were 
generic, it would imply a relation between the m remaining values, p 2D ~ m ^ . . . , P, and the 
previous values, R 2 , . . . ,p 2 - D-m+1 ; since an arbitrary state is not subject to such a restriction, 
the POVM could not be PS I-complete. We conclude that if a measurement is to be PS I- 
complete, a transversal intersection at every step in our process is generically necessary. In 
accordance with Sard's theorem, a nontransversal intersection at any step could be removed 
by slightly perturbing the values of our constraints. 

Now consider the final step in our process. As we have just argued, if the POVM is to be 
PS I-complete, this step generically involves the transversal intersection of a compact, one- 
dimensional submanifold, X 2 d-2, with the hyperellipsoid Sp, which is defined by Eq. (J2UJ) 
(with P > generically), and the result of this final intersection is a set of points, 2 2 d-i- 
The submanifold X2D-2, being compact and one-dimensional, is a set of closed curves that 
do not intersect one another or themselves. Moreover, the hyperellipsoid Sp is closed and 
orientable, thus splitting M 2D_1 into an inside and an outside. A closed curve intersecting 
such a surface transversally must intersect it an even number of times, since any entrance 
from the outside must be paired with an exit. Thus we can conclude that I2D-1 consists 
generically of an even number of points. 

We need, of course, a stronger result than this because the inversion symmetry, £ — * — £, 
already implies that the points in X2D-1 come in pairs. What we need to show is that 
I2D-1 contains at least two points with £0 > 0, but we can get this by a simple extension 
of the above reasoning, which is illustrated in Fig. |21 Let X® be the part of the hyperplane 
£o = that lies inside or on Sp. The inversion symmetry guarantees that the points in the 
intersection of I2D-2 with X Q come in pairs, £ and The points in a pair are distinct, 
because an intersection at £ = is forbidden by the normalization constraint. 

Now consider the union of the £ > part of Sp, which we denote as S p + \ and X . Since 
this union is closed and orientable, we can conclude, just as we did above for Sp, that the 
closed curves in I2D-2 intersect it an even number of times. Since we have just established 
that I 2 d-2 intersects X an even number of times, we find that T 2 d-2 intersects S p + ^ an 
even number of times. This means that the final set, I2D-1, generically has an even number 
of points with £ > 0, corresponding to distinct pure states consistent with the outcome 
probabilities. This completes the proof of Theorem 1, because this generic ambiguity, at 
least two-fold, shows that a POVM with 2D — 1 elements cannot be PS I-complete. 

IV. SUFFICIENCY OF 2D OUTCOMES 

In this section we construct an explicit example of a PS I-complete POVM with 2D 
outcomes, thus showing the sufficiency of this number of outcomes for PS I-completeness. 
Combined with the result of the preceding section, this shows that a minimal PS I-complete 
POVM has 2D elements. 
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Letting {\j), j = 0, . . . , D — 1} denote an orthonormal basis for a D-dimensional quantum 
system, we define the operators 

X jk EE \j){k\ + \k){j\, 

Y jk ee -i|j)<*l+i|*>0'l, (21) 

EE |j><j|-|*X*l- 

These operators can be thought of as the Pauli operators in the subspace spanned by \j) 
and | A;}. In this section we only need the Oj versions of these operators, but we use others 
in Sec. El 

We now consider a POVM consisting of the following positive operators: 
a|0)<0| , 

b(l + X oj ) and b(I + Y 0j ), j = 1, . . . , D — 1, (22) 
T . 

Here a and b are positive numbers, and T is a "throw-away" POVM element that must 
be included to satisfy the completeness condition (JTJ. We can always make T a positive 
operator by choosing a and b small enough. Notice that this POVM has 2D elements, as 
promised. We now demonstrate that this POVM is PS I-complete. 
We let 

D-l 

M>> = I>|j> (23) 

3=0 

be an arbitrary normalized pure state, expanded in the basis We remove the global 

phase freedom by choosing c = r to be real and nonnegative, and we write Cj = Xj + %y 3 - 
for the remaining amplitudes. 

The first outcome probability for the POVM (}2*2*|) is 



p° = a\(0\ij)\ 2 = arl , (24) 
from which we get the (positive) amplitude of the state |0), 

ro = M ■ (25) 
V a 

The remaining outcome probabilities (except for the throw-away outcome) are 

p x * = + X oj )m = b(l + 2r oXj ) , 

P" J = Hm + YojM = b(l + 2r oyj ) , (26) 

for j = 1, . . . , D — 1. Except for states with r = (a set of measure zero), these, too, can 
be immediately inverted to give the remaining amplitudes, 

yj 2br n K ' 
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This completes the (trivial) demonstration that the POVM (}2*2*|) is PS I-complete. 

The inversion procedure fails on the (D — 1) -dimensional subspace orthogonal to |0), 
where r = 0. There is a simple way to handle this failure in a tomographic procedure, 
i.e., in a sequence of measurements of this POVM, made on many identical copies of the 
same system, with the outcome probabilities estimated from the outcome frequencies. All 
one has to do is to precede the tomographic procedure with a single premeasurement in 
the basis Whatever the outcome of this premeasurement, that outcome has nonzero 

probability, so we can choose it to be the state |0). 

Notice that after determining r , the reconstruction of the remaining amplitudes Cj is 
reduced to solving linear equations. This transformation of a fundamentally quadratic in- 
version into a linear one is the reason it is so easy to show that the POVM (}2*2*j) is PS I- 
complete. The key to this transformation is that the probabilities f!26|) have no terms that are 
quadratic in the amplitudes Cj. This, in turn, is a consequence of having the unit operator in 
the "middle" POVM elements, b(I+ X j) and b(I+Y j), since for normalized states, these unit 
operators always put a 1 into the corresponding outcome probabilities. We can't just leave 
the unit operator out, because we have to put something there to make the middle POVM 
elements positive. Suppose, for example, that we tried to make the middle POVM elements 
rank one by changing them to b(Poj + X j) and b(Poj + Y j), where Poj = |0)(0| + is the 
projection operator onto the subspace spanned by |0) and \j). Then the probabilities (|2T)j) 
become p x ^ = b(r^ + x? + y? + 2r Q Xj) and p y, i = b(r^ + x 2 - + + 2r yj). The quadratic 
term x| + introduces a generic two-fold ambiguity into the determination of each of the 
amplitudes Cj, which makes the POVM not PS I-complete. 

We should stress that the role of the POVM (j22j) is solely to demonstrate that 2D 
outcomes are sufficient for PS I-completeness, because the presence of the unit operator 
in the middle POVM elements makes this POVM a poor candidate indeed for an actual 
tomographic procedure. The unit operators mean that the middle outcomes all have nearly 
the same probability, with only a weak dependence on the amplitudes Cj. It is from this weak 
dependence that one must extract the amplitudes, and this makes a tomographic procedure 
based on this POVM very inefficient. 

Weigert [10] has described a quite different way of making a PS I-complete measurement, 
which we can manipulate to produce a single PS I-complete POVM with 2D outcomes. 
Considering the system to be a spin-s particle (D = 2s + 1), Weigert first shows that 
the probabilities for two spin components, S ■ n and S ■ n', along infinitesimally different 
directions, are sufficient to determine a generic pure state up to 2 2s = 2 D ~ 1 ambiguities. He 
then shows that the ambiguities can be resolved by knowing just the expectation value of 
the spin component S ■ e along the direction e orthogonal to n and n'. 

To describe this scheme, let \n, m) denote the eigenvector of S ■ n with eigenvalue m. 
Notice that 



is a potential POVM element whose expectation value contains the same information as 
the expectation value of S • e. Weigert describes getting the expectation value of S ■ e in 
the standard way from measurements of that spin component, but we can get the same 
information directly by including F in a POVM. Furthermore, since the measurements of 
S n and S n' are separately normalized, we can amalgamate two projectors, one from each 
measurement, into a single POVM element without losing any information. If we put all this 
together, Weigert's work shows that a generic pure state is determined by the expectation 



F = — 



(5 ■ e + si) 




(29) 



m=—s 
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values of the operators 

\n,m)(n,m\ , m — — s + 1, . . . , s, 
\n',m)(n',m\ , m = — s + 1, . . . , s, 
\n, — s)(n, —s\ + \n', —s){n', —s\ , 

F . (30) 

This list contains 2D POVM elements, whose expectation values are PS I-complete, but 
the list does not constitute a POVM because the elements don't sum to I. To turn the list 
into a POVM, we use a trick mentioned in the Introduction in the context of IC-POVMs. 
Suppose we have a set of positive operators, {F c }, which are PS I-complete in that their 
expectation values are sufficient to reconstruct a generic unnormalized vector up to 
a global phase. First form the operator G = ^2 C F C . This operator is invertible, for if it 
weren't, there would be a vector |e) such that (e|G|e) = 0, which would imply that F c \e) = 
for all values of c ; this would mean that the expectation values of the operators F c provide 
no information about the amplitude (e^), so the operators would not be informationally 
complete. Now define new operators E c = G~ X I 2 F C G~ X I 2 . These operators are positive and 
sum to I, so they make up a POVM; moreover, they are obviously PS I-complete. It is also 
worth noting that if F c is rank-one, so is E c . 

Although it is easy to calculate G for the list (j3*Uj) . it is not very enlightening to calculate 
the POVM elements. We simply note that the trick can be applied to the list ()30|) . thereby 
turning Weigert's PS I-complete measurements into a minimal PS I-complete POVM — and, 
moreover, one with 2D — 2 rank-one POVM elements. The price of this transformation is 
that the measurement no longer corresponds to just measuring a few spin components. 

One might be tempted to apply this trick to the POVM elements (|22|) . without the throw- 
away T, since the probability for outcome T never appears in the inversion procedure. We 
know this can't work, however, because the result would be a PS I-complete POVM with 
2D — 1 elements, contradicting the results of Sec. IIHI The reason it doesn't work is that the 
inversion procedure assumes normalized pure states, and the trick only works if the inversion 
works for normalized and unnormalized states. For unnormalized states, the throw-away 
can't be discarded because it is required to fix the normalization. 

V. RANK-ONE PS I-COMPLETE MEASUREMENT WITH 3D - 2 OUTCOMES 

In this section we consider PS I-complete POVMs whose POVM elements are all of rank 
one, i.e., multiples of one- dimensional projectors. We don't know the minimal number of 
elements in a rank-one PS I-complete POVM. We suspect it is close to or even equal to 
2D — the reworked Weigert measurement discussed in Sec. IIVI has 2D elements, with just 
two not of rank one — but the best we can do for the present is 3D — 2 elements, an example 
of which we present now. 

The setting here is the same as in Sec. IIVI with states and operators defined as in Eqs. (J22J) 
and (J21|) . except that we now allow the states to be unnormalized. Consider the following 
four states in the subspace spanned by \j) and \k), 

IVW = \j) , 

fei) = cos(0/2)|j)+sin(0/2)|A:) , 
Ite) = cos(0/2)|j)+e 2 "/ 3 sin(0/2)|fc), 
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fe 3 ) = cos(0/2)|j) + e- 2m /^m{e/2)\k) , (31) 
where we assume < 9 < it. Notice that we can write 

IVW = cos(0/2)|j> + e 2 ^ a - 1 )/ 3 sin(^/2)|A;) , a = 1,2,3. (32) 
The corresponding one- dimensional projection operators are given by 

Zjk) i 



sin ^ + y ifc J^sin^ + cos 



|fe;o)(feo| 




Pjk 


1^1) (fell 




Pjk 


|^jfc;2)(^jfc;2| 










>- 



; -sin0 - F ife 



sin 6 1 + Zjjn cos 6> J . (33) 



When cos# = —1/3, these are the four tetrahedral states in the two-dimensional subspace 
spanned by \j) and \k); in this subspace, they are specified by the Bloch vectors (JTTJ). Notice 
that 



X>iU(feal = ^(Pjk + Z jk cos6) =3(cos 2 (6/2)\j)(j\ + sm 2 (6/2)\k)(kty . (34) 

a=l 

We now consider the following POVM elements, numbering 3D — 2, 
a|0)(0|, 

b\il>j-ij. a )(il>j-ijia\ , j = 1, • • • , -D — 1, a = 1,2,3, (35) 

where we assume a, b > 0. The first thing to show is that these POVM elements are sufficient 
to reconstruct a generic unnormalized pure state \ip). We begin by determining cq = r > 
from the probability for the first POVM element, a|(0| , ?/')| 2 = clTq. As we work through the 
remaining amplitudes, suppose that we have determined c , . . . , Cj_i and are now trying to 
get Cj from the three POVM probabilities 

&|(^--i,i;a|^)| 2 = %os(fl/2) Cj _i + e- 27ri(Q - 1)/3 sin(#/2) Cj | 2 , a = 1,2,3. (36) 

Writing these three probabilities out, it becomes clear that they determine Cj, provided 
Cj_i 7^ 0. The only hitch in this method is that we get stuck if any Cj except the last is 
zero, so we have a set of D — 1 (D — l)-dimensional subspaces where the method fails. This 
failure set is a set of measure zero, however, so the procedure succeeds in reconstructing a 
generic pure state. 

We now use the trick that we applied to Weigert's measurement in Sec. I1VI We form the 
sum of the operators in Eq. (|35j). 

D-l 3 

G = a|0)(0| + X)X) 6 l^-^«><^-^«l 

j=l a =i 

D-2 

= (a + 36cos 2 (^/2))|0)(0|+36^|j)(j|+3&sin 2 (^/2)| J D-l)(D-l| , (37) 

3=1 



15 



where we use the sum (|34|) . Choosing 6=1/3 and a = sin 2 (6'/2), we get 

G = I-cos 2 (0/2)|L>-l)(L>-l| . (38) 

The trick introduced in Sec. IIVI now tells us that the following POVM elements make up 
a POVM that is PS I-complete: 

sin 2 (#/2)G- 1/2 |0)(0|G- 1/2 = sin 2 (fl/2)|0>(0| , 

^G-^l^-ijjaXVi-ijjalG- 1 / 2 = l\ij>j-ij;«)(ipj-ij;«\ , 3 = h ■ ■ • , D - 2, a = 1, 2, 3, 
^G-^IVd-^-i^)^^^-!^^- 1 / 2 , a = 1,2,3. (39) 

Here 

G-^I^D-lja) = CO S (0/2)|£> - 2) + e ^(«-l)/3|£) _ 1) , a = 1, 2, 3. (40) 

It is easy to see directly, without appealing to the general validity of our trick, that these 
particular elements make up a rank-one POVM that is PS I-complete. 

Notice that D = 2 is a special case, because the sum in the second form of Eq. ()37|) is 
missing. By making the tetrahedral choice of angle with a = b = 1/2, we get the four- 
element tetrahedral POVM of Eq. ([lip. Thus the POVM (|3*9*jl can be considered as a sort 
of generalization of the tetrahedral measurement. 

Amiet and Weigert [13] have formulated a PS I-complete measurement procedure for a 
spin-s particle, which is based on measuring the probabilities for spin components along any 
three noncoplanar axes. By amalgamating three of the outcomes, one from each of the spin 
components, this procedure can be reduced to a PS I-complete POVM with 3D — 2 elements, 
but with the one amalgamated element not of rank one. 

VI. CONCLUSION 

To determine any pure state of a D-dimensional quantum system, except perhaps a set of 
measure zero, requires a measurement with at least 2D outcomes. That is the key result of 
this paper. We also demonstrated that a 3D — 2 element POVM composed only of rank-one 
operators (multiples of one- dimensional projectors) is sufficient for determining a generic 
pure state. A number of open questions suggest themselves. Is it possible to construct a 
2D-element POVM that can determine any pure state, without having to exclude sets of 
measure zero? If not, how many extra elements must be added to achieve this stronger form 
of PS I-completeness? What is the minimal number of rank-one POVMs necessary for PS I- 
completeness? How do the results obtained here connect with PS I-completeness for wave 
functions of infinite-dimensional quantum systems? Finally, we have done next to nothing 
in this paper on the question of the efficiency of PS I-complete POVMs for determining pure 
states from outcome frequencies in an actual tomographic procedure. We know that the 2D- 
element PS I-complete POVM introduced in Sec. IIVI is very inefficient, but we have nothing 
to say at present about what efficiencies can be achieved by minimal PS I-complete POVMs. 
These questions suggest avenues for research that might turn PS I-complete POVMs into 
useful tomographic tools in situations where identical systems are confidently thought to be 
in a pure state. 
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VII. EPILOGUE 

2D, or not 2D: That was the question. Nature, her pride to assuage, tempts with 
another: whether 'tis nobler in the mind to suffer the slings and arrows of 3-D — 2 rank-one 
POVM elements or to take arms against a sea of equations, and by opposing, solve them? 
Alas, Nature hath her secrets hidden well; yea, even the mightiest among us hath failed her 
veil to penetrate. Be not discouraged!, dear reader, for we leave to thee the Herculean task 
to cast out these D — 2 elements, banishing them thus, as so many angels fallen from the 
Kingdom of Heaven. States without States, measuring without Measurements. Parvo non 
ex nihilo [T3|. And if, upon completion, 2D carries the day: Lo! such beauty to behold! 
Ne'er shall our POVMs a finer form take, for it be optimal. And should 3D — 2 remain, 
recall: "A dream itself is but a shadow." A shadow of what? "... a shadow's shadow." 
Thus the Symphony of the Universe plays: movement upon movement, endlessly subtle, till 
time itself reigns alone as King. 
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